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ANALYSIS OF THE THREE-POINT-BEND TEST FOR MATERIALS WITH 
UNEQUAL TENSION AND COMPRESSION PROPERTIES 
by Christos C. Chamis 
Lewis Research Center 

\ 

SUMMARY 

An analysis capability is described for the three -point -bend test applicable to ma- 
terials of linear but unequal tensile and compressive stress-strain relations. The capa- 
bility provides numerous equations of simple form and their graphical representation. 

The various equations are derived using linear structural mechanics principles. Pro- 
cedures are described for determining the local stress concentration in the vicinity of 
the load point and the failure mode at the critical stress location. 

Examples are given to illustrate the usefulness and ease of applicability. The func- 
tional behavior of the equations is described graphically. This is done to expedite the 
analysis and to illustrate the influence of material properties and geometric parameters 
on properties such as: maximum deflection, maximum tensile (compressive) stress, 
maximum shear stress, flexural modulus, shear modulus, and local stress concentra- 
tions. 

Comparisons are made with values obtained when the material is considered to have 
equal tensile and compressive mechanical properties. Some typical results from these 
comparisons follow: The maximum bending deflection and stress can be underestimated 
by 25 percent or more. The flexural modulus can be underestimated by 35 percent, and 
the shear modulus by 25 percent. The horizontal shear stress can be underestimated by 
15 percent in the vicinity of local stress concentrations. No corrections are required for 
the maximum bending stresses and deflection when the length -to -depth ratio is 20 or 
greater. Failure initiates on the tensile side when the material’s compressive strength 
is equal to or greater than its tensile strength. 

The capability can also be used to reduce test data from three -point -bend tests, ex- 
tract material properties useful in design from these test data, select test specimen di- 
mensions, and size structural members. 



INTRODUCTION 


Structural resins are used extensively as matrices in advanced fiber composites. 
These resins have moduli and strengths that are different in tension and compression. 
Typical ranges of ratios are, for the compression -modulus to tensile -modulus ratio, 

1. 0 to 1. 5 and, for the compressive-strength to tensile strength ratio, 1. 5 to 2. 0. 

Structural resins are generally characterized by means of the three-point -bend test. 
The test data are reduced using simple beam formulas that do not account for different 
tensile and compressive properties. The reduced data can be in considerable error 
depending on the magnitude of the difference between the tensile and compressive prop- 
erties. 

There are several good reasons why the three -point -bend test is used extensively in 
material characterization: economy, simplicity of specimen preparation and testing, 
ease of adaptability to enviromental testing, suitability for cyclic loading and fatigue 
testing, convenience for fracture toughness studies, and the availability of well docu- 
mented simple formulas for analyzing materials having equal tension and compression 
properties (refs. 1 to 4). Another important reason is that the three -point -bend test is 
a simple way to subject a specimen to tension, compression, and shear simultaneously. 
In this sense, a three -point -bend test provides a direct measure of the structural in- 
tegrity of the material. 

Available simple beam formulas are used for analysis because no comparable 
formulas are available which account for different tensile and compressive mechanical 
properties. Therefore, an investigation was performed to derive all the equations 
needed for the analysis and test data reduction of three -point -bend test for materials 
with unequal tensile and compressive mechanical properties. The specific objectives 
were the 

(1) Derivation of all the equations needed in the form of simple formulas 

(2) Investigation of the stress concentration effects in the vicinity of the load point 

(3) Examination of failure stress and failure initiation. 

The governing equations were derived using well known linear structural mech- 
anics principles. Symbols are defined when they first appear and are also compiled in 
appendix A for your convenience. The equations are described in the main text where 
several examples are given to illustrate their simplicity and usefulness. The detailed 
derivations are presented in appendix B. Although the terms ’'bending” and ’’flexural" 
are used interchangeably, herein, I use only flexural for consistency. Analogous equa- 
tions can be derived for the four-point -bend test by following the procedure used in 
appendix B. 
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THEORY AND USEFUL EQUATIONS 


Theoretical Background 

Closed form equations for analyzing or reducing data from a three -point -bend test 
(fig. 1) are described herein. The equations are for shift of the neutral plane (zero 
bending stress plane), maximumum bending deflection, flexural (bending) modulus, 
shear modulus, maximum bending stress, and maximum shear stress. Only the final 
equations are given herein. The detailed derivations are given in appendix B. 

The equations are derived following the procedure used in classical beam theory but 
accounting for different moduli in tension and compression. Briefly, classical beam 
theory is based on the elastomechanics concepts of force equilibrium, stress-strain re- 
lations, strain displacement relations, and the conservation of elastic strain energy. In 
the present development, it is assumed that both the stress-strain and strain displace- 
ment relations are linear. 

The following equations and relations were used to develop the theory presented 
herein (refer to fig. 1). The force equilibrium equation at a point (x, z) in the beam is 
given by 


^£-^ 1=0 ( 1 ) 
dx dz 

where a and r denote normal and shear stress, respectively. The resultant force 
at a section x is given by 



( 2 ) 


where h^ is measured from the neutral plane to the bottom of the beam and h^. to the 
top (fig. 1(b)). The subscripts T and C denote tension and compression, respectively. 
The resultant moment at a section x is given by 


M = b 



(3) 


where b is the beam width. The linear stress strain relations are given by 

= E»p^ijt = T = Gy (4) 
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where e and y denote normal and shear strain, respectively, and E and G denote 
normal and shear modulus, respectively. 

The linear strain displacement relation at a point z is given by 



(5) 


2 2 

where w is the bending deflection and (d /dx ) denotes the second derivative. 

The details on how equations (1) to (5) were used to derive the governing equations 
of interest herein are given in appendix B. The property and governing equations are 
briefly described in subsequent sections. The governing equations are used to illustrate 
the effect of the modular ratio E T /E C on the property under discussion. 


Neutral Plane Location and Bending Stiffness 
The location of the neutral plane (plane of zero bending stress, fig. 1(b)) is given by 



( 6 ) 



( 7 ) 


where h is the beam thickness (fig. 1). 

The shift of neutral plane as a function of modular ratio E t /E c is presented 
graphically in figure 2. As can be seen in this figure, the neutral plane shift exceeds 
10 percent when the modular ratio is less than 0. 8 or greater than 1. 2. It is useful to 
note from equation (7) and figure 2 that the neutral plane shifts towards the higher- 
modulus material. 

The bending stiffness, also known as bending rigidity or flexural rigidity, is given by 

D = 3 ( h b E T + h t E c) ( 8 ) 
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where D denotes the bending stiffness. The effects of the modular ratio on the bending 
stiffness will be illustrated later. 

Equation (8) can also be given as 


or 


or 


D-IbhhX 


(8a) 


D = 


A bhh?E 
3 1 


(8b) 



(8c) 


It can be readily verified by inspection that equations (8) reduce to the well known ex- 

O 

pression Ebh /12 when the modular ratio equals one. 

Equation (8c) indicates that when the tensile modulus only is used to calculate bend- 
ing deflections, the deflection will be underestimated if E T /E ( -, <1.0 and overestimated 
if E^E^ > 1. 0. 


Maximum Flexural Deflection 

The equation for computing the maximum bending deflection (deflection at midpoint) 
is given by 


Pi 3 

w = — 
max 48D 

where P is the load, l is the length, and G is the shear modulus. 

The shear modulus for a material with different tensile and compressive moduli is 
given by (see appendix B) 


|1 + 1. 6 bhpi 



(9) 
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( 10 ) 



where tVp is the Poisson’s ratio measured on a tensile test. 

The influence of the modular ratio on the flexural deflection is shown in figure 3 
where the variable w max / w p raax is plotted against l / h for various Er^/E^ . The 
variable w jr max is the flexural deflection for the case where E^/E^, =1.0 and for no 
shear contribution. It is given by 


w Fmax 



4bh 3 E T 


( 11 ) 


Two points are worthy of note in figure 3: 

(1) The modular ratio has significant influence on the flexural deflection. The 
simple formula (eq. (11)) underestimates the flexural deflection by 25 percent when 
E t /E c =1.5 and overestimates it by 25 percent when Erp/E^ = 0. 5. These observa- 
tions were anticipated from the discussion following equation (8c). 

(2) The length -to -depth ratio 1 / h has negligible effect on the flexural deflection 
when 1 / h > 10 in the range of modular ratio 0.5< E^/Eq ^ 2. 0. This is an important 
observation because the researcher does not have to be overly concerned about the 
specimen length so long as it has a length -to -depth ratio greater than 10. 


Bending and Shear Moduli 


The apparent flexural modulus from a three -point -bend test is obtained by solving 
equation (11) for E. The result is 


E F = 


P V 


4bh 3 w 


( 12 ) 


Fmax 


where E p is the apparent flexural modulus and w Fmax is the measured deflection 
under the load. Equation (12) assumes that the neutral plane coincides with the mid- 
plane of the beam. 
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The error introduced in equation (12) when Erp/E^, ^1.0 may be calculated by 
using w Fmax from equation (9) in equation (12). This is illustrated in figure 4 where 
the nondime nsional variable Ep/E T is plotted against the modular ratio for i / h = 10 
and 40 (flexural modulus from deflection curves). As can be observed from figure 4, 
the influence of the modular ratio on the apparent flexural modulus is considerable. 

The flexural modulus may also be estimated by the approximate relation, which is 
derived by averaging the stiffnesses (see appendix B) 



where Ep denotes the approximate flexural modulus. The ratio Ep/E T is plotted 
against E T /E ( ~, in the figure 4 curve labeled ’’flexural modulus” (approximate). As 
can be seen in figure 4, the approximate flexural modulus is a good approximation (less 
than 5 percent error) when 0. 75 < E^/Eq < 2. 0. 

The effects of the modular ratio on the shear modulus are obtained from equa- 
tion (10). The shear modulus from equation (10) nondimensionalized with respect Gjqq 
is plotted against E^/E^ in figure 5, where Gjqq is the shear modulus of an isotropic 
material and is given by 


G ISO~ 


E 

2(1 + v) 


(14) 


As can be seen in figure 5, the modular ratio E^/E^ has considerable effect (as high 
as 25 percent) on the shear modulus. 

The conclusions to be made from the previous discussion and the curves in figures 
4 and 5 are 

(1) Modular ratio E^/E^, has significant influence on the flexural and shear moduli 
(could be as high as 35 and 25 percent, respectively, for 0. 5 < E^/E^ < 2. 0) 

(2) The length -to -depth ratio i / h has insignificant effects on the flexural modulus 

(3) The flexural modulus may be approximated sufficiently closely using a simple 
equation. 


Maximum Flexural and Shear Stresses 

The maximum flexural tensile stress occurs at x= i/ 2 and z = -h^ (fig. 1(b)). 
The equation to predict this stress is given by 
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a Tmax 



(15) 


or 


Tmax 


3P _/l) 

4bh b ^h y 


(15a) 


The maximum bending compressive stress occurs at x = l/2 and z = h^. (fig. 1(b)). 
The equation to predict this stress is given by 


a Cmax 


fc 

4D 


(16) 


or 


Cmax 


-3P 1 1 
4bh t \h / 


(16a) 


Letting h^ = h/2 in equation (15a) and h^. = h/2 in equation (16a) result in 


a 

max 


1. 5 


EL 

bh 2 


(17) 


which is the equation given for a material with equal moduli in tension and compression. 
Nondimensionalizing equation (15a) with respect to equation (17) results in 


a Tmax h 
a max ^b 


(18) 


The corresponding result for 0Q max is given by 


a Cmax 

a max 



(19) 


Equations (18) and (19) indicate that the maximum flexural stress (tensile or compres- 
sive) occurs on the extreme fibers of the stiffer side. The graphical representation of 
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°Tmax and CT Cmax a S ainst E T /E C is shown in figure 6, which shows that the cor- 
rections to the maximum stresses can be as much as 20 percent. 

Two more useful stress ratios are obtained by dividing equation (15a) by equa- 
tion (16a) and then using equation (6). The results are 


g Tmax _ M 1/2 

CT Cmax h b \ E C / 


( 20 ) 


Equation (20) indicates that the ratio of the maximum stresses varies as the square root 
of the ratio of their respective moduli. 

The maximum shear stress is given by the equation 


T xzmax 



( 21 ) 


which is the same as that given for a material with Erp/E^ =1.0. The reason for this 
condition is that the shear stress equals zero at the top and bottom surfaces of the beam. 
These boundary conditions can only be satisfied by a parabolic variation of the shear 
stress through the beam thickness. However, the plane of maximum shear stress is at 
z = 0 and depends on the modular ratio. This dependence is illustrated graphically in 
figure 2. 

The important points to be noted from the previous discussion are 

(1) Simple relations for various stress ratios exist 

(2) The stress ratio a T max /°c max varies as the square root of the modular ratio 
and attains its largest magnitude on the extreme fibers of the stiffer side 

(3) The correction to either ^ Tmax or tf Cmax can be as large as 25 percent when 
compared with a material with equal moduli 

(4) The magnitude of the maximum shear stress is independent of the modular ratio. 
However, the plane at which the maximum shear stress occurs depends on the modular 
ratio. 


STRESSES AND DISPLACEMENTS IN THE VICINITY OF THE LOAD 

The equations described previously predict the stresses quite accurately at a dis- 
tance one-half the thickness away from the load application point. In the vicinity of the 
load these equations need additional corrections. In this vicinity (point A, fig. 1(b) both 
x and z direction stresses have the same sign, and equations with equal moduli are 
applicable. 
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In reference 5 approximate equations are presented for a material with equal 
moduli for tension and compression. The equations from reference 5 are presented 
herein and discussed briefly to illustrate the effects of stress concentration. 

The stress in the x-direction in the section under the load is given by (see fig. 

7(d)) 






where 


ll < z < — 


and for z = -(h/2) 



( 22 ) 


(23) 


(24) 


Both equations (22) and (24) indicate that the stress concentration correction to ct x 
is insignificant and becomes negligible as if h > 10. A graphical representation of the 
local stress concentration contributions to a x through the beam thickness is shown in 
figure 7(a). Note that the correction decreases both the tensile and compressive stress 
values predicted by the beam formula. Stated differently, the beam formula slightly 
overpredicts the flexural stresses. 

The stress in the z-direction (fig. 7(d)) at the section under the load is given by 


a z = 


nbh 


1 + 2 



(25) 


where 




( 26 ) 
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and 


as z — - 
2 


(27) 


<7 


Z 


Note that the condition indicated by equation (27) is typical of elastic solutions near 
singularities. However, a_ will produce permanent damage locally when a has ex- 
ceeded the compressive strength of the material. 

The graphical representation of equation (25) is shown in figure 7(b), where the cr 
variation through the thickness is plotted. As can be seen in figure 7(c), a decays 
very rapidly away from the load point. 

The shear stress variation through the beam thickness at a section in the vicinity of 
the load (fig. 7(d)) is shown in figure 7(c). This figure was generated from data given in 
reference 5 (p. 105). The interrupted portion of this curve indicates the lack of data in 
that region. The maximum shear stress occurs at z = h/4 (fig. 7(c)) with an approxi- 
mate value of 


t V7 « 0. 86 1- (28) 

xz bh 

The maximum shear stress from the simple beam formula is given by equation (21). 

The ratio of equation (28) to equation (21) is 

?HJ ocal = 0^86 = h 15 (29) 

t 0 75 

xz simple -beam 

This leads to the conclusion that the maximum shear stress can be about 15 percent 
greater than that obtained by means of the simple beam theory and that it occurs at about 
z = h/4. 

The maximum deflection is also affected by the stress concentration in the vicinity 
of the load. For a material with equal moduli in both tension and compression, the 
equation from reference 5 (p. 107) is 
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For a material with equal moduli, in both tension and compression, equation (9) re- 
duces to 


w = 


pr 


4bh 3 E 


1 + 1.2 




(31) 


Assuming v= 0. 4 and using the isotropic relation E = 2(1 + v) G, the middle term in 

2 

equation (30) becomes 3. 00(h/z) . The corresponding term in equation (31) is 
o 

3.36(h/i) . The conclusion from this exercise is that equation (31) and, therefore, 
equation (9), overestimates the shear contribution to the maximum deflection. It can be 
verified by direct substitution that both shear and local stress contribute less than 
1 percent to the maximum deflection when (i/ h) is 20 or greater. 

The previous discussion leads to the following conclusions: 

(1) The local stress concentration has negligible effect on the bending stresses and 
bending deflections predicted by simple beam formulas when l / h is equal to or greater 
than 20. 

(2) The through-the -thickness normal stress o_ is very high in the immediate 
vicinity of the load point. It decreases very rapidly as the distance away from the load 
point increases. 

(3) The local stress concentration contributes about 15 percent to the maximum 
shear stress predicted by the simple beam formula, and it occurs at about the top one- 
fourth point. 


EFFECTS OF STRESS CONCENTRATION ON THE SPECIMEN FAILURE STRESS 


The effects of the local stress concentration on the specimen failure stress are con- 
sidered in this section. The stress state at point A (fig. 1(a)) is biaxial compression. 

An equation describing the combined-stress-state failure envelope for a material with 
different strength in tension and compression is given by (ref. 6 with K^ 2 an< * K ^ 12o?/3 
replaced by K XZ( ^) 


F ( a x’ ff z’ a xz) = 1 - 





■ K XZQf)3 


g XQ? g zj3 
S xa S x/3 


(32) 


where F is the combined stress failure function, a is the stress due to load, S is the 
unaxial failure strength, K depends on the elastic constants of the material and is unity 
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for isotropic materials, the subscripts x and z denote direction, and the subscripts 
a and /3 denote tension or compression. 

The envelope described by equation (32) is shown in figure 8 for a material with 
room -temperature properties representative of epoxy resins and a value of K- Xza p = 1. 0 
Combined or uniaxial stress states causing failure are represented by points on the 
failure envelope curve in figure 8. 

Applying this information to the failure stress of a three-point -bend test and recall- 
ing the local stress concentration effects on the flexural stresses from the previous sec- 
tion lead to the following conclusion: Failure will initiate on the tension side of a three - 
point -bend test specimen that is made from materials with compressive strength equal 
to or greater than the tensile strength. 

The preceding conclusion is obvious from both intuition and figure 8 for the case 
when the compressive strength is greater than the tensile. The case of equal strengths 
is a little more subtle and can be shown by the following procedure: 

(1) Assume a representative value of if h. Herein we assume if . h = 20. 

(2) Assume the specimen fails on the tension side (point B, fig. 1(a)) when 

°xB “ ®xT' 

(3) Calculate the load P required to produce a stress <x x g = S x ,p. 

(4) Knowing this load, calculate the corresponding stresses o and a in the 
immediate vicinity of point A. 

(5) Plot this point in the tension -tension quadrant in figure 8. The tension quadrant 
is sufficient since we are interested in the case of equal compressive and tensile 
strengths; that is, S xC " S xT- 

The point plotted in step (5) should fall inside the region bounded by the failure envelope 
to verify the conclusion for the case of equal strengths. 

Following the preceding procedure yields: From equation (24) and ifh. = 20, 


From the simple beam formula cr 

A 


Equation (25) at 2z/h = 0. 95 yields cr z bh/p « 25.4. Using P from equation (33), we 
obtain 


p fts r!l s 
30 


'xT 


(33) 


= (1.5 Vlfb\T) (from fig. 8(a)) and P (from eq. (33)): 


°xA “ °- 957 S xT 


ff zA= °- 


85 S. 


'xT 
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The combined stress and the single stress are plotted as points A and B in figure 8. 
Point A in figure 8 is inside the failure envelope as was previously postulated. There- 
fore, the conclusion is verified for the equal strength case. 

The l/h at which a zA = a xA is of interest. This ratio may be obtained readily 
following the previous procedure. The required equations are 


1.5 P l . 1.30 ,, 
= + * 


xA 


(34) 


bh 


bh 


ff zA 



Equating these two equations yields 



(35) 


(36) 


The interpretation of equation (36) is as follows: 


> 17 ®xA increases faster than a ^ 



a xA increases a t the same rate as <t z a 



a xA increases slower than 


It is noted that as - l/h approaches one, cr zA will increase considerably faster than 
and produce local bearing failure under the load. At l/ h= 1, a zA /cr xA ~ ® f rom equa- 
tions (34) and (35). 

The important conclusions from the previous discussion follow: 

(1) Failure initiates on the tension side of a three -point -bend test specimen. This is 
the case for specimens made from materials with compressive strength equal to or 
greater than the tensile strength and with an l/h greater than 17. 

(2) The apparent flexural strength will be equal to or greater than the smaller of the 
simple strengths. If S xT < S xC , then S xT < flexural strength < S xC - 
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ADDITIONAL APPLICATIONS OF THE EQUATIONS AND GRAPHS 

In addition to analyzing the three -point -bend test, the equations and graphs pre- 
sented herein can be used in other ways. Some examples are briefly described. If the 
tensile modulus and the flexural modulus are known, the corresponding compressive 
modulus can be obtained from figure 4. The shear modulus can be obtained from fig- 
ure 5. The neutral plane location is obtained from figure 2. If, in addition, the appar- 
ent flexural failure stress is known, the actual flexural failure stress may be obtained 
from figure 6. 

The corresponding simple tensile strength and a lower bound on the compressive 
strength may be estimated by the use of equations (15a) and 16(a) and figure 7(a). 

If a specimen fails in compression and the apparent flexural strength is known, the 
actual flexural strength can be obtained from equation 16(a). The corresponding simple 
compressive strength may be estimated using equation (16a), figures 7(a) and (b), and 
equation (32). 

The actual and apparent flexural strengths can be estimated when the simple tensile 
or compressive strengths are known. The estimate may be obtained by following the 
procedures used in the previous section. 

Equations (34) and (35) can be used to specify a loading rate ratio of 
(fig. 1). H the specimen is suspected to have failed by a combined stress state, includ- 
ing shear, then this condition can be checked by computing the stresses cr , a , 
at that point and substituting in equation (32). 

It is worth noting at this point that the procedures described herein apply to fiber 
composites as well. However, the data in the figures do not because the shear modulus 
was determined from equation (10) or (14), neither of which is valid for fiber compo- 
sites. The procedures described previously are also applicable in sizing structural 
components. 


SOME COMMENTS ON INSTRUMENTING THE THREE -POINT -BEND TEST 

The equations locating the neutral plane may be used in conjunction with figure 7(c) 
to locate strain gages for measuring the shear modulus using a three -point -bend test. 

A direct determination of the location of the neutral plane can be obtained by mea- 
surement. This is accomplished by placing strain gages on the top and bottom surfaces 
of the specimen (near points A and B, fig. 1) along the x-direction. The gages should 
be placed at a distance h/2 away from the load to avoid local stress concentration 
effects. 
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Plotting the stresses predicted by equations (15a) and (16a) against the strains of the 
corresponding gages yields the required stress-strain curves which can be used to de- 
termine both the tensile and compressive moduli. 


SUMMARY OF RESULTS 

Results obtained from analyses of three -point -bend tests for materials with different 
moduli in tension and compression are summarized in this section. The tensile- 
modulus to compressive -modulus ratio is referred to herein as the modular ratio (MR). 

1. An analysis capability was generated. The capability can be used to analyze the 
test data or select test specimen geometry for three -point -bend tests. 

2. The maximum flexural deflection is sensitive to MR. The simple beam formula 
underestimates the maximum deflection by 25 percent when the MR = 1. 5 and over- 
estimates it by 25 percent when MR = 0. 5. The maximum bending deflection is insen- 
sitive to the length -to -depth ratio if h when l/h is greater than 10. 

3. The MR has significant influence on the flexural and shear moduli. Corrections 
of 35 and 25 percent, respectively, are needed for the simple beam formula predictions 
when MR = 0. 5. But l/h has negligible influence on the flexural modulus. 

4. The flexural stress ratio (tensile to compressive stress) varies as the square 
root of MR. The maximum stress magnitude in tension (compression) occurs in the ex- 
treme fibers of the stiff er side. The correction to maximum tensile (compressive) 
stress can be as large as 25 percent compared to a material with equal moduli. 

5. The local stress concentration in the vicinity of the load has negligible influence 
on the maximum flexural stress and deflection where l /h equals 20 or greater. How- 
ever, the maximum shear stress is increased by about 15 percent and occurs at about 
the one -fourth depth point from the top of the beam. 

6. Failure will initiate on the tension side of a three -point -bend test specimen made 
from material with compressive strength equal to or greater than the tensile strength. 

7. The equations and graphs can be readily used to extract data which are useful for 
design purposes and as a guide to instrumenting the test specimen. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 8, 1973, 

501-22. 
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APPENDIX A 


SYMBOLS 


b specimen, beam width 

D flexural stiffness 

E normal modulus 

F combined stress failure function 

G shear modulus 

h specimen thickness 

h^ distance from neutral plane to 

bottom 

h t distance from neutral plane to 

top 

K xza)3 coef f icient in combined stress 
failure function 

l specimen length 

M resultant moment 

P applied T .oad 

S simple strength 

w bending deflection 

x coordinate 

z coordinate 

y shear strain 


e normal strain 

v Poisson's ratio 

o normal stress 

r shear stress 

Subscripts: 

b bottom 

C compression 

F property obtained by simple 

beam formulas 

ISO isotropic material property 

max maximum value 

T tension 

t top 

x direction along x-axis 

z direction along z-axis 

a T or C, tension or compres- 

sion 

/3 T .or C, tension or compres- 

sion 
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APPENDIX B 


DERIVATIONS 


Basic Relations 


The linear uniaxial stress-strain relations for a material with different moduli in 
tension and compression are given by 

0" rp — 

(Bl) 

°c = e C E C 

(B2) 

r = y G 

(B3) 


The strain-displacement relations of interest in these derivations are given by (see 
fig. 1(d)) 


e = -z 


d 2 w 

dx 2 


(B4) 


The force equilibrium equations of interest in these derivations are (figs. 1(a) and (b)) 


do 3r 
x _ xz 

9x 3z 


= 0 


(B5) 


where 0 < x < f and h^ < z < h^, 


f° f h t 

L a T dz + J a - 0 


(B6) 


where 0 ^ x < i , and 
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(B7) 


M = b 



zarpdz + b 



za c dz 


where 0 < x < l . 


Equation of the Neutral Plane 


Using equation (B4) in equations (Bl) and (B2), substituting the result in equa- 
tion (B6), and carrying out the integration yield 



f„ 


X t ! o d2w 1 9 

o^dz = - E m h“ ^ - - ErK 


>T“b ^2 


^ ^2 


= 0 


' From which follows 


and 



(B8) 


(B9) 


Equation (B9) is the equation of the neutral plane. It can also be expressed in terms 
of h. U sing the identity 


h = h, + h b 


(BIO) 


in equation (B9) and collecting terms yield 



(Bll) 
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Equation of Flexural Rigidity (Stiffness) 

Using equation (B4) in equations (Bl) and (B2), substituting the result in equa- 

2 2 

tion (B7), and carrying out the integration noting that d w/dx is independent of z 
yield 



h ‘W 

0/ dx 2 


From which follows 


\ d\_ 

/dx 2 

2 2 

Let the coefficient of d w/dx be D. Thus, 

D -|( h b E T + ,1 t E c) 


M 


“ -J b (^ E T 


W E C 


(B12) 


(B13) 


Equation (B13) is the flexural rigidity of a material with different moduli in tension and 
compression. Equation (B13) can be expressed in terms of either h^ or h t and h 
through the use of equation (B8) or (B9) and equation (BIO): 


D = J bhh b E T = J bhh t E C 


The resultant moment in terms of D is 


(B14) 


M= -D 


d 2 w 


dx 


(B15) 


Stress Moment Relations 
Multiplying by z both sides of equation (B15) yields 
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Mz = -Dz — 
dx 2 


and using equation (B4) results in 

Mz = De 


Substituting equations (Bl) and (B2) in the last equation and rearranging yield, re- 
spectively, 


Em 

CTm = MZ — 
1 D 

0 < x < i -h b < z < 0 

(B16) 

and 



E c 

c D 

-4-* 

VI 

tsj 

VI 

o 

(SJ 

VI 

X 

V 

o 

(B17) 

But in the region 0 < x < (l/2) 

1 

II 

s 

(B18) 

Using equation (B18) in equations (B16) and (B17) yields, respectively, 


„ Pxz _ 

Um — Em 

T 2D T 

0 < x < I -h, < z < 0 

2 b 

(B19) 

»c= E c 

u 2D U 

0 < x < i- 0 < z < h, 

2 1 

(B20) 

The maximum tensile and compressive stresses are obtained by the following substi- 
tutions: 

x = 1 and 
2 

z = -h b in equation (B19) 


x = - and 
2 

z = h t in equation (B20) 
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The results are, respectively, 


ff Tmax = E T 


and 


P£ h t 

CT Cmax “ 4D E C 

Using the values of D from equation (B14), the last two equations become 

= JLL_ L 

Tmax 4bh ^ h 

_ = I 

Cmax 4bh ^ h 


(B21) 


(B22) 


Letting h b = h t = (h/2) in equations (B21) and/or (B22) results in 


■^Cmax a Tmax 


3PZ 

2h 2 


(B23) 


which is the well known maximum stress formula. Equations (B21) and (B22) constitute 
the maximum stress formulas of interest in this investigation. 


Shear Stress Applied Load Relation 
From equation (B5) we have 




(B5) 


In the region 0 < x < (i/2) and 0 < z < h t , cr x is given by equation (B20). Substituting 
equation (B20) in equation (B5), carrying out the differentiation with respect to x, and 
integrating the result with respect to z yield 
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Pz* 

4D 


Applying the boundary condition at z = h^ 


-t xz + f (x) = 0 

and r _ = 0 yields 
xz 


Ph t 

F(x) = — - E r 
4D c 

Substituting this value of F(x) in the previous equation and rearranging yield 



(B24) 


The maximum shear stress is obtained by letting z = 0 in equation (B24), using the 
value of D from equation (B14), and using equation (BIO). The result is 


T xzmax 


3P 

4bh 


(B25) 


Equation (B25) is the same as that of a material with equal moduli. However, the maxi- 
mum value is not at the midplane of the specimen. 

Carrying out the same procedure for the region 0sx<(j/2) and -h^ < z < 0 
yields 


T xz = 


ph b E t 

4D 


1 - 



(B26) 


But at z = 0 the shear stress predicted by equation (B26) must equal that of equa- 
tion (B24). Carrying out the algebra and reducing the result yields 
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The last equation defines the plane of maximum shear stress. Note that this equation 
is the same as equation (B9), which is the equation for the neutral plane. Therefore, 
the neutral plane and plane of maximum shear stress coincide as is the case for ma- 
terials with equal moduli. 


Maximum Flexural Deflection 


The total energy stored in the beam is given by 



Using the stress-strain relations (eqs. (Bl) to (B3)) and the stress-force relations 
(eqs. (B19), (B20), and (B26)) in equation (B27) and carrying out the integrations and 
simplifying yield 


U = 


l b 


48D 


30GD 


4 h b + E C h f; 


(B28) 


The midpoint deflection (maximum deflection) is obtained by differentiating equa- 
tion (B28) with respect to P (Castigleano’s first theorem). The result after simplifica- 
tion is 


w 


max 


-w - p * 

' W '/ 2 ' 48D 


1 + 


1.6bh2hE 2 /h b 


GD 




When E t = Eq, h^ = h^ = (h/2), equation (B29) reduces to 



(B29) 


(B30) 


Equation (B30) is the simple beam formula when the shear contribution is included. 
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The shear contribution to the maximum deflection may be limited to specified values 
by adjusting the U/h) ratio in equation (B30). For example, if the shear contribution is 
to be limited to 1 percent or less* then l / h is 

> ^120E^ 1/2 

For isotropic materials E/G « 3 and 1 / h > 19. However, for advanced fiber com- 
posites for which equation (B30) applies, E/G « 30 and i / h > 60, which is considerable. 
This procedure can be used to size components from these materials. 



Relations Between Shear, Tensile, and Compressive Moduli 

Relations for G and E,p and E^-, are obtained in the following manner: Refer to 
figure 9(a). For the stress state shown, the strain energy stored in a unit -volume - 
element is 



T ct T 


C C T 


: C U C 


Using equations (Bl) and (B2) yields 


U = 



E, 


a C°T 



(B31) 


in terms of stresses only. 

For the stress state in figure 9(b) the energy stored in a unit -volume-element in 
terms of stresses only is 


U = 



V T 

+ —°T°C 
E T 



(B32) 


Assuming no energy loss during the stressing -straining process, the energy is the same 
in both cases. Subtracting equation (B31) from equation (B32), we obtain 
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(B33) 




If in figures 9(a) or (b) we let the magnitude of = ct t , we have a case of pure shear 
that can be expressed as 

t= |^£i| = tfip (B34) 

The energy stored in a unit -volume -element for the case of pure shear is given by 



2G 


2 


T 


(B35) 


But the energy in equation (B35) is equal to that of equation (B31) or equation (B32) since 
it is produced by the equivalent final stress states. Equating equation (B35) to, say, 
equation (B31), we obtain 


2G 


T 2 = 


2 Erl 




E, 


' °C°T + 


2Er 


(B36) 


Using equations (B33) and (B34) in equation (B36), we obtain, after simplification, 



(B37) 


Equation (B37) reduces to the well known case for E^ = E T as can be verified by in- 
spection. 

When E t , PVp, and G are known from independent tests and not from the relation 
E = 2(1 + v)G, then the compressive modulus E^-. can be computed from equation (B37). 
The result is 



G 


E 


T 


- (1 + 2i/ t ) 


26 



Equation (B38) indicates that a material with Erp = G(1 + 2i/rp) is incompressible. 

We note here, that analogous derivations leading to relation of equation (B37) are 
presented in reference 7. 


Approximate Equation for the Flexural Modulus 

Assume that the flexural modulus may be approximated by the average extensional 
modulus through the beam thickness. In equation form it is 


E 


F 


“ r (h b E T - W 




(B39) 


Using the relation E ( -,/E t 



(eq. (B8)) in equation (B39) and simplifying yields 


E 


F 



(B40) 


as an approximate equation for determining the flexural modulus. 
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(a) Specimen geometry. 


Compress- 



lb) Normal stress (c) Shear stress Id) Normal strain 
variation. variation. variation. 


Figure 1. - Schematic of geometry and stress variations of 
three-point -bending test specimen from material with un- 
equal tensile and compressive properties. 



Neutral 

plane 


Modular ratio, E-p/E^ 


Figure ^ - Effect of modular ratio on shift of neutral plane 
for three-point bending. 
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Normalized flexural modulus, Ep/Ey . / q \/_,q 

r 1 Normalized maximum deflection, w^frEyj/Pr 


Modular ratio. 

e t ,e c 



Figure 3. - Maximum deflection of three-point-bending beam 
made from material with unequal tensile and compression. 
Moduli maximum deflection normalized with respect to that 
of equal moduli. 



Figure 4. - Effect of modular ratio of flexural modulus for three- 
point bending. 




Modular ratio, Ey/E c 


Figure 5. - Shear modulus of material having 
different values for tension and compression 
moduli; Poisson's ratio - 0.4. 


P 



Figure 6. - Effect of modular ratio on maximum bending 
stress for three-point bending a /b ■ 40). 
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(c) Peak shear 
stresses at 
0 < x < h/2. 


id) Schematic near vicinity of 
load application. 


Figure 7. - Stress variation in vicinity of applied load for modular 
ratio, Ej/E c , 1. (For geometry notation part (d) this figure. ) 
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Stress, ksi 

Figure 8. - Failure envelope for typical epoxy resin in 
biaxial stress field. 
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<b) Compression then tension. 

Figure 9. - Stress application for determining relations between 
shear, tensile, and compressive moduli. 
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